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The proton's gluon structure 
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The proton's gluon structure function at small x is larger than nowadays is commonly believed. 



1. Introduction 

We do not understand perturbative QCD at 
small x. In particular, as I will explain, we 
do not understand how to apply DGLAP evo- 
lution there. However, if we combine it with 
Regge theory and use an important message from 
the HERA data for the charm structure function 
F£(x,Q 2 ), it is possible|l[| reliably to extract the 
gluon structure function g(x, Q 2 ) at small x. It 
turns out to be larger than nowadays is commonly 
believed. This is seen in figure 1. The most recent 
CTEQ and MRST structure functions @|| agree 
well with each other and with those extracted by 
the two HERA experiments 0] because they all 
use similar procedures; however, Donnachie and 
I believe that the old MRSG structure function is 
nearer the truth. 

2. Regge theory — the two pomerons 

When one tries to fit data, it is usually sensible 
to start with the simplest assumptions and then 
refine them later. In its simplest form, Regge 
theory leads to fixed powers of x at small x, and 
it turns out that two terms are enough: 



F 2 (x,Q 2 ) ~ h{Q 2 )x- C0 + h(Q 2 )x~ 



(1) 



The second term corresponds to soft-pomeron ex- 
change, with ei w 0.08 determined from soft re- 
actions. The data need a term that rises more 
rapidly at small x; one needs eo ~ 0.4. By fitting 
the data at each Q 2 , Donnachie and I found j|] 
that a successful and economical parametrisation 
of the coefficent functions is provided by 

fo(Q 2 ) = A (Q 2 ) 1+e °/(l + Q 2 /Qo) 1+£o/2 
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/i(Q 2 )=Ai(Q a ) 1+ «/(l + Q a /Q?) 1+ei 



(2) 



with Q » 3 GeV and Qx w 0.8 GeV. To make 
the fit, we used real-photon data and DIS data 
with x < 0.001, so that Q 2 ranges from 0.045 to 
35 GeV 2 . If we then simply multiply the resulting 
form (|I]) by (1 — x) 7 , as is suggested by the di- 
mensional counting rules HQ , it agrees quite well 
with the HERA data even beyond x = 0.1 and 

rto Q 2 = 5000 GeV 2 . This is shown in figure 
Note that this factor (1 — x) 7 should not be 
taken too seriously; it is much too simple. 

Data[|| for the charm structure function 
F£(x, Q 2 ) have the remarkable property that, at 
all available Q 2 , they fit to just the single hard- 
pomcron power of a;. Further, to an excellent ap- 
proximation the coupling of the hard pomeron 
appears to be flavour blind: 



F^(x,Q 2 ) = f c (Q 2 )x~ 



with 



f c (Q 2 ) = 0.4 f (Q 2 ) 



(3) 



(4) 



So if we define a charm-production cross section 



a c (W) 



47r 2 a EM 



F 2 c (x,Q 2 ) 



x=Q 2 /(W 2 +Q 2 ) 



(5) 



it behaves as W 2e ° at all Q 2 , even down to 
Q 2 = 0: see figure ||. Perturbative QCD directly 
relates F% {x, Q 2 ) to the gluon structure function, 
so that at small x it too must be dominated by 
hard-pomeron exchange alone, even at quite small 
values of Q 2 . This what causes the rapid rise at 
small x of the DL curve in figure [l]. 



2 



35 
30 
25 

15 
10 

5 





0.0001 



x DL 






02 8 


MRSG"\ 




CTEQ6M 




MRST2001 



0.001 



Figure 1. Gluon structure functions 

from http : / /cpt 1 9 . dur . ac . uk/hepdat a/ pdf 3 . html 
and reference 1 
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Figure 3. Data|| for the electroproduction of 
charm at various Q 2 , with W°' 8S and pQCD fits 
(upper and lower curves, respectively) 
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Figure 2. Regge fit to ZEUS and HI data for 
F 2 (x,Q 2 ) for Q 2 between 0.045 and 5000 GeV 2 . 
The parameters were fixed using only data for 
x < 0.001 and therefore Q 2 < 35 GeV 2 . 



Figure 4. NLO evolution of the hard-pomeron 
coefficient function (solid curve) with the phc- 
nomenological fit (broken curve) 
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3. DGLAP evolution 

The singlet DGLAP equation introduces the 
two-component quantity 

where t — log(Q 2 /A 2 ). If we Mellin transform 
with respect to x, the equation becomes very sim- 
ple: 

-u(N,Q 2 )=P(N,a s (Q 2 ))u(N,Q 2 ) (7) 

The usual approach is to expand the splitting ma- 
trix P in powers of a s . However, this is mathe- 
matically illegal when N is small. Compare 

y/N 2 + a s -N = a s /2N - a 2 s /8N 3 + . . . (8) 

Here, each term in the expansion is singular at 
N = but the function itself is regular there: 
the expansion is illegal for iV 2 < a s . Similarly, it 
is likely that whenever expanding P(N, a s (Q 2 j) 
makes it large, and therefore makes u(N, Q 2 ) vary 
rapidly with Q 2 , the expansion is dangerous. 

My own belief is that P(N, a s (Q 2 )) has no sin- 
gularities in the complex- TV plane, or at least no 
relevant singularities. My reason is that solving 
(0) would cause a singularity of P(N,a s (Q 2 )) to 
induce an essential singularity in u(N, Q 2 ) (that 
is, a nasty one). The variable N is closely re- 
lated to the orbital angular momentum I, and 
I was brought upjj| to believe that matrix el- 
ements such as u(N, Q 2 ) do not have essential 
singularities in the complex Z-planc. This point 
of view contrasts with that of those who believe 
that the value of eo is associated with a singular- 
ity of P(N, a s (Q 2 )) and may even be calculated, 
perhaps by refining the BFKL approach. I think 
that very probably eo is a nonperturbative quan- 
tity that therefore cannot be calculated. 

A fixed-power behaviour x~ e ° of F^{x, Q 2 ), 
such as in (|l|), corresponds to an A^-plane pole: 

< 9 » 

If we insert this into the DGLAP equation (]?]) 
and equate the residue of the pole on each side of 
the equation, we find 

^f(Q 2 ) =P(N = e ,a s (Q 2 ))i(Q 2 ) (10) 



eo is far enough from for the expansion of 
P(N = e ,a s (Q 2 )) to be reasonably safe. So we 
may easily use the DGLAP equation to calculate 
the evolution of the hard-pomeron component of 
F2(x,Q 2 ). But this is not the case for the soft- 
pomeron component, because e± « 0.08 is too 
close to 0. 

According to figure [l], the various gluon struc- 
ture functions come together at x rj 0.01. It is 
reasonable to assume that for values of x larger 
than this the evolution of the two elements of 
u(x, Q 2 ) does not use values of N close to 
and therefore the conventional analysis is correct. 
So we can start at some not-too-large value of 
Q 2 , 20 GeV 2 say. We determine the value of 
fo(Q 2 ) there from the phenomenological fit (j|) 
and f g (Q 2 ) from the MRST gluon structure func- 
tion xg(x,Q 2 ), which for x greater than about 
0.01 fits very well to x~ £ °(l — x) 5 . We choose 
Anlo such that a s {M'j c ) = 0.116 and use ( |l0| ) 
to calculate^] the NLO evolution of fo(Q 2 ) an d 
fg(Q 2 ) m both directions. The result for fo(Q 2 ) 
is the continuous curve in figure |. The dashed 
curve is the phenomenological form (J2[). Provided 
we adjust A so that still a s (M^-) = 0.116, LO evo- 
lution gives almost identical results; this is shown 
in figure ^|. 

The agreement between the pQCD calculation 
and the phenomenological curve is a success not 
only for the concept of the hard pomeron, but also 
for pQCD itself. The evolution is from a single 
value of Q 2 , not the customary global fit[^,||], and 
it introduces far fewer parameters. 

Notice that, as Q 2 increases the large- a; be- 
haviour of xg(x,Q 2 ) becomes steadily steeper 
than (1 — a;) 5 , and so the largest value of x for 
which x~ e ° is a good approximation to the struc- 
ture function steadily decreases. Figure ^ shows 
an estimate of this. 

We may use the gluon structure function to cal- 
culate the charm structure function F^x, Q 2 ). 
The result, using just LO photon-gluon fusion 
with a charm-quark mass m c = 1.3 GeV, is the 
solid curves in figure [?| This is an important 
check on the consistency of the approach. As 
is seen in figure 0, a steep gluon distribution is 
needed to fit the data at small Q 2 . 
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In conclusion, the conventional approach to 
evolution needs modifying at small x. It can 
be corrected if we combine it with Regge theory, 
but only partly — we can only treat the hard- 
pomeron part. The resulting gluon distribution 
is larger at small x than has so far been supposed 
and gives a good description of charm production. 
I should add that we want good data for the longi- 
tudinal structure function, because this gives the 
most direct window on the gluon distribution. 
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Figure 7. Charm structure functions at Q 2 — 1.8 
GeV 2 , with ZEUS datag 
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Figure 5. LO and NLO evolution of the hard-pomeron and gluon coefficient functions 
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Figure 6. Gluon structure function at Q 2 — 20 and 200 GeV 2 



